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Recently single molecule experiments have shown the importance of internal fric¬ 
tion in biopolymer dynamics. Such studies also suggested that the internal friction 
although independent of solvent viscosity has strong dependence on denaturant con¬ 
centration. Recent simulations also support such propositions by pointing out weak 
interactions to be the origin of internal friction in proteins. Here we made an attempt 
to investigate how a single polymer chain with internal friction undergoes reconfigu¬ 
ration and looping dynamics in a confining potential which accounts for the presence 
of the denaturant, by using recently proposed “Compacted Rouse with internal fric¬ 
tion (CRIF)”. We also incorporated the effect of hydrodynamics by extending this 
further to “Compacted Zimm with internal friction (CZIF)”. All the calculations 
are carried out within the Wilemski Fixmann (WF) framework. By changing the 
strength of the confinement we mimicked chains with different degrees of compact¬ 
ness at different denaturant concentrations. While compared with experiments our 
results are found to be in good agreement. 

I. INTRODUCTION 

In recent past there have been efforts based on single molecule experiments to elucidate 
the role of internal friction in protein folding BE]. These single molecule experiments showed 
internal friction to play a significant role in folding especially when the protein starts with a 
more compact unfolded state. Among the different types of experiments carried out to study 
the dynamics of the unfolded proteins, the most common one is a combination of FRET and 
nsFCS, where two residues n and m of a protein are labelled with fluorescence probes and the 
fluctuation of the distance between them is observed from the efficiency of energy transfer. 
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An auto-correlation function is then calculated from the distribution, which eventually gives 
a characteristic relaxation time, also called the reconhguration time (r„m) [T]. In another 
type of fluorescence experiment, two different residues of a protein are again tagged with a 
fluorescence quencher and donor. The time required for the the donor to be quenched which 
happens only within a certain distance, gives an estimation of loop formation time {Tnm,ioop) 
in a protein [3]. Both the time scales, Tnm and Tnm,ioop seem to have a nonzero intercept when 
plotted against the solvent viscosity. This residual reconhguration or looping time has been 
then attributed to “internal friction” which is present within the protein and assumed to be 
independent of the solvent around. Unfortunately a molecular picture of this internal friction 
is still lacking but it is believed that dihedral rotations, hydrogen bonding and other weak 
interactions contribute to internal friction. Recent simulations on proteins support such 
propositions as well gHZ]. This also tells why internal friction is more when the protein is in 
more compact unfolded state. Though the notion of “internal friction” is not very old in the 
physical chemistry community it has always been a topic of research among polymer 

rheologists. For example Rabin and Ottinger looked at the origin of internal viscosity in 
a Gaussian chain mm- Following an idea of de Gennes na they derived an expression 
for the relaxation time, Trei associated with internal viscosity as Trei = /kBT{r]s -|- r]i) 
where, R = aN'^ and a, N are the monomer size and chain the length respectively, v is the 
Flory exponent [IMS]. Therefore in the limit solvent viscosity r]s —)■ 0, it has a non-zero 
intercept proportional to the internal viscosity rji. This is in the same spirit as that of 
Khatri and McLeish [TB], where a modihed Rouse model gives a mode relaxation time which 
is dependent on internal friction, -|- Tint- Such a model gives a reconhguration 

time between any two monomers n and m as Tnm — 0.82r^ -|- Tint UTl, where t^ is the 
rouse time, proportional to the solvent viscosity. Although such a model qualitatively can 
predict the trend of the variation of reconhguration time as a function of solvent viscosity 
and produce a non-zero intercept but cannot account for the changes in “internal friction” at 
diherent degrees of compactness encountered in experiments. It is expected that the protein 
experiences diherent level of “internal friction” at diherent concentrations of the denaturant 
as the compactness of it changes. This can be seen from the plots of the reconhguration time 
[2] against the solvent viscosity at diherent denaturant concentrations where the intercepts 
gives the time scale due to “internal friction”. The higher the denaturant concentration the 
smaller the intercept. In the present study we extend recently proposed models [niiiB] to 
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include the effect of the compactness of the polymer chain to internal friction by introducing 
a confining harmonic potential to each monomer of the chain which accounts for the change 
in chain conformation due to denaturant. At lower denaturant concentration chain is more 
compact so the confining potential is steeper but as the denaturant concentration increases 
the confining potential becomes shallower. With this model, the looping dynamics is studied 
within Wilemski Fixman(WF) framework [19] assuming the polymer chain to be Gaussian. 
Loop formation between any two parts of a bio-polymer is supposedly the primary step of 
protein folding, DNA cyclization [201 El]- It is worth mentioning that WF formalism seems 
to work fine here and has been used to calculate the same in presence of hydrodynamic 
interactions by Chakrabarti and the effect of viscoelastic solvent [231121]. The method 
has also been used in the past extensively [251 - 130] to investigate other aspects of the polymer 
looping problems. 

The paper is arranged as follows. The details of the the polymer models are given in 
section II. Section III deals with the methods. Results and discussions are represented in 
section IV. Section V is the conclusion. 

II. POLYMER MODEL 

In the Rouse model, a polymer chain is imagined as a series of brownian particles con¬ 
nected by harmonic springs with no hydrodynamic interactions and excluded volume effect 
[3IIE2]. If Rn{t) is the position of the monomer at time t, where n can be varied from 
0 to V, then dynamics of Rouse chain with [N + 1) monomers is described by the following 
equation of motion 


dRn{t) ^ 

^ dt 

Where C, is the solvent viscosity and k = 
f{n,t) is the random force with moments 


^ +/W>G (1) 

is the spring constant with kuhn length b. 


ifiriR)) = 0, {fain,ti)ffiimR 2 )) = 2C,kBT5c,fib{n - m)5{ti - G) 


( 2 ) 


Standard procedure to treat such a system is to describe the dynamics in the form of 
normal modes {Xp) 
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^j^dXp{t) _ 


dt 


= -K^pi^) + fpi^) 


(3) 

rR 


Here, = 2N( for p > 0 and (q = N(. The relaxation time of mode is = ^ = 

F F f^p i' 

and kp = ^^ 2 ^ ■ The slowest relaxation time = 3 ^ 2 called Rouse time. 

Rouse model does not take into account the effects of hydrodynamic interactions. The 
simplest possible model which takes care of it, is due to Zimm [311 Eg. It is possible to show 
that in 6 condition under pre-averaged hydrodynamic interaction, the Zimm chain can be 
described by the same Eq. ([^ but with a different scaling of friction co-efficient with the 
mode number Cp = C ^ where, kp remains the same i.e. k^ = kp = kp^^"> [22]. The 


relaxation time of mode is r^f = ^ = 


p3/2 


and = 




Q^^3/2j^^T ’ slowest relaxation 

time or Zimm time. 

Recently in polymer rheology community internal friction has been introduced within 
Rouse description using a dashpot between the neighbouring monomers. This is commonly 
referred to as Rouse with internal friction (RIF) which has the following euation of motion 


dRn{t) / d\ d'^Rnit) 


(4) 

Khatri and McLeish [T6| showed that above model can also be treated with normal mode 
description as is done with Rouse chain. Similarly a Zimm chain with internal friction 
(ZIF) can be constructed. In both the models RIF and ZIF, internal friction appears as an 
additive constant to the relaxation times of each normal modes (p > 0) [T71 [33|. The new 
relaxation time is = Tp + Tint where. Tint = Qnt/k. This is because of the redehned 


qR(z)if ^ + 27r‘^p‘^(int/N,p > 0 where 


remains 


unchanged, i.e. = k^^^\ 


"p '"p 

In RIF and ZIF there is no scope to introduce the effect of denaturant concentration to 
the chain conformation. It is obvious that the more is the denaturant concentration, the 
less compact the chain is and the less is the internal friction. To incorporate such effect, one 
possibility is to put all the monomers in a conhning potential. This becomes analytically 
trackable when the conhnement is harmonic. Thus the restoring force acting on the 
monomer would be ~ 0)^) where, kc is the spring constant. Similar models of 

polymers in conhned potential has been used in other contexts too [341139] . This model is 
named as Compacted Rouse with internal friction (CRIF). We would also like to mention 
that essentially the same model has recently been used to £t the simulated results [g. 
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No Hydrodynamics 


Parameter 

Rouse 

RIF 

CRIF 

Cp 

Cp^ = 2iVC 

Cp^^ = Cp+ 2F'^p‘^Cint/N 

^CRIF^^R^2FVCint/N 

kp 

pR _ Gir’^ksTp^ 

“ FP 

uRIF _ u 

hjp — /ip 

= 2Nkc + kp 

Tp 

— Tp 

P kp p2 


^CRIF _ Cp 1 ^CRIF 

T “ kCRiF -r ^int 

'^int 

0 

= Cint/k 

^CRIF _ tint 

int k-\-kcN^ 


TABLE I: List of parameters for Rouse, RIF and CRIF 




(5) 


Similarly ZIF can be extended to Compacted Zimm with internal friction (CZIF). Now, 
as well as are mode dependent becanse = kp^^^ + 2Nkc, bnt the 

friction coefficients associated with normal modes remain nnchanged i.e. cind 

qCzif _ (^ziF^ ^ conseqnence, the relaxation time for mode of this Ronse chain 


-ZIF 

bp Sp 

is, 


W^F + fc_|_fc ^ 2 /p 2^2 and for this Zimm chain in 9 solvent nnder pre-averaged 
hydrodynamic interactions is Tp^^^ = • The time scale for the internal 

friction of CRIF as well as CZIF become identical Depending on 

the valne of kc, mode nnmber p and the chain length two extreme sitnations can arise. 
For example if kc is small, chain is short (small N) and p » 1, k + kcN"^/ p^f'^ can be 

Same is trne for Zimm chain. This is 
expected as this corresponds to a very weak conhnement of the chain where higher modes 
rarely contribnte. On the other hand, if k^ is large, then k + kcN"^ /p^f‘^ — k(.N‘^ /p^f‘^ and 
Therefore in snch a sitnation k^. controls the dynamics of the polymer 
and higher modes contribnte more. All the parameters for Rouse and Zimm chain are 
depicted in Table P and [IT| respectively. 


approximated as k and then 
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Hydrodynamics 


Parameter 

Zimm 

ZIF 

CZIF 

Cp 


Cf^ = C/ + 27rVCWA 

C^CZIF ^^z 

kp 

> ~ FP 

l-ZIF _ u 

tvp — rvp 

= 2Nkc + kp 

Tp 

II 

II 

= ip + 

^CZIF _ c# , CZIF 

T “ kCziF -T T^nt 

'^int 

0 

Tin/ = Cint/k 

^CZIF _ Cm 

int k+kcN'^ 


TABLE II: List of parameters for Zimm, ZIF and CZIF 

III. CALCULATION METHODS 

A. Reconfiguration time and mean square displacement 

The time correlation function of normal modes for Rouse as well as Zimm chain is given 
by the following equation [31] 

7 rjn 

{Xpa{0)Xqfj{t)) = -^Spq6al3ew{-t/rp) ( 6 ) 

r£p 

One needs to only choose the correct Tp to find the correlation for Rouse or RIF and Zimm 
or ZIF [T8] . 

When normal mode correlations of CRIF or CZIF are constructed the basic structure of 
correlation functions remain the same. Only kp and Tp need to be replaced by and 

tCR{z)if respectively. 

(Ap„(0)Xg^(f)) = exp (7) 

Kp 

Therefore Eq. ([^ is a generalized expression. From now on, all the expressions are defined 
for compacted polymers with internal friction to make the discussion more extensive, from 
which one can easily get back to RIF and ZIF by taking the limit kc —t 0. Usual Rouse and 
Zimm models can also be recovered by putting fcc = 0 and Qnt = 0 

OO 

The distance between and monomers is Rnm{t) = 2 X) — cos(2^)], 

p=i 

OO 

where R„(t) = Aq + 2 X) Ap(t)cos(^) and has the following time correlation function 

p=i 
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0nm(t) = {Rnm{^)-Rnm{t)) = 4 £ [cOS ( ^ ) - COs{^^)f exp{-t / (8) 

Reconfiguration time Tnm is calculated by taking an time integration of (f>nm{t) [IZl HH SO]. 

OO 

T~nm J ( 9 ) 

0 

Where, (j)nm{t) is the normalized correlation function which is defined as 




4^nm (^) 


In the limit when is high, Tnm = N 


2 p=i 


( 0 ) 

C + C*„tp2 7r^/]V2 




E ^Jcos{E^)-cos{I^)]^ 

P =1 

dent of k as expected and dynamics is controlled by kc- 

Similarly, MSD of R„m(^) has the following closed form expression. 


( 10 ) 

, which is indepen- 


OO 

{Rnm{t) - -Rnm(O))^^ = 4 ^ 


6knT 


, CR(Z)IF 
p=i kp 


cost 


,p7rn^ 

If' 


— COS[ 


^pirm 

If 


( 11 ) 


B. Looping time 

A well known procedure to calculate the looping time is due to Wilemski-Fixmann (WF) 
[T^ . Although the method was initially posed for end-to-end loop formation, but can be 
generalized for looping between any two monomers of the polymer chain [THl 123] , as long as 
the chain is Gaussian. So, this method is applicable to CRIF and CZIF as well. Looping 
time within WF framework is given by following expression, 

r„„.,„p= (12) 

Jo / 

Where, Cmnf) is the sink-sink correlation function between the and monomer 


Cnmf) = J dR nm J dRyim^ oS{R 

nm )G{R 

nm? t\R 

nm,0? 0)S{R 

nm,0 )PiR 

nm,0 ) (13) 












G{Rnm)t\Rnmfi)^) IS the Conditional probability that the distance between and 


monomers of the chain, is Rnm,o at time t = 0, and Rnm at time t. The chain is assumed to be 


3/2 




in equilibrium at t = 0 with distribution, P{Rnm,o) 


exp 


S{Rnm) is the sink function having inhnite strength which depends only in the separation 
between the two monomers which basically takes care of the loop formation between two 
monomers. This delta function sink is used in many cases to calculate the survival probability 
of a brownian particle IHHH]. 

The conditional probability is given by. 



3/2 


(14) 


Choosing the sink function S{Rmn) to be a delta function the expression for looping time 
becomes 



exp[-2xo^l^{i)/{l - ^l^{t))]sinh[{2xo^nrn{t))/{I - 0L(^))] 
(2Xo0nm (t) ) Jl- 


(15) 


Where, 


?>a? 


(16) 


^ \^^mn/eq 

For CRIF and CZIF (R^o)e(j following analytical expression 



(17) 


In the limit, kc —)■ 0, 



(18) 


Obviously {Rnq)^^ > {RNo)eq conhning potential makes the polymer more compact 

and as a consequence the effective kuhn length becomes shorter. Thus the equilibrium 
distribution function of CRIF and CZIF have narrower widths but remain Gaussian [2]. 
Similar situation of renormalization of the kuhn length arises during swelling of a polymer 
chain [45], where the kuhn length gets longer. 
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IV. RESULTS AND DISCUSSIONS 


Schuler and his group studied the dynamics of cold shock protein (Csp) using FRET 
and FCS. From the auto-correlation function of the distance separating any two arbitrary 
residues they determined intramolecular diffusion coefficients [a 06]. In our analysis all 
the parameters have been chosen in accordance with these experiments. Thus the polymer 
has 67 monomers {N -|- 1 = 67), with kuhn length b = 3.8 x 10“^°. In consistance with 
the viscosity of water, the solvent friction ( = 9.42 x at temperature 30077. 

The force constant of the springs connecting monomers to each other is calculated from the 
relation k = SksT/b"^. Without a clear molecular picture of the internal friction, we have 
to invoke an ansatz which is Cmt = Cint,of{nb) [IE], where nb is the number non-adjacent 
monomers contributing to the internal friction and f{nb) = (cq + CiHb -|- C 2 uf -|- ....), where 
(int,o is the internal friction due to adjacent neighbour interactions as shown for RIF in Fig. 
Q which arises if Cq = 1 and Cj = 0 for i > 0, Ub = 0. As the solvent quality drops, the 
polymer shrinks bringing more and more number of monomers together. Hence, nb and c* 
increases so as f{nb). Taking this relation one step further, we associate k^ with Qnt with 
another ansatz kc = kc + -|- ciUb + C 2 U? -|- ....) where, A is typically in the order of 

TT^/V^. Now kc has two parts kc and the second part which is proportional to (^int- This is 
obvious since the more compact the polymer is, the higher is the internal friction. Again in 
the limit c* = 0 for i > 0 and Uft = 0 /Cc,o = k^ accounts for the situation when 

the chain experiences a harmonic conhnement but does not have internal friction. We have 
considered kc to be zero in our calculations for the sake of simplicity. 

and identical and referred to as in general. Therefore in Fig. 

@ the time scale due to the internal friction of each normal mode or is plotted 

against normal mode number p at three different values of kc namely kcp, 27^,0) 3/Cc,o- As 
expected it saturates at higher modes which shows that higher modes get rarely affected by 
the conhnement and to the higher modes adds as an additive constant as mentioned 

earlier. It can also be seen from the graph is completely independent of the mode 

numbers when kc is zero since in the absence of kc, is nothing but . With 

increasing kc, (int also increases which is why becomes larger. 

The generalized expression for the reconhguration time between two monomers m and 
n is given in Eq. (g, with m = 0 and n = N it becomes the end-to-end reconhguration 
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time. End-to-end reconfiguration time is computed at different solvent viscosities [r]) with 
three different values of for CRIF. Effects of hydrodynamics have also been taken care 
of in CZIF. To compare with experiments, reconfiguration times are plotted against rj/rjo 
in Fig. (|^ and Fig. Q where, rj is solvent viscosity and rjo is viscosity of water. In 
experiments solvent viscosity is generally varied by adding viscogenic agents such as glycerol 
to the solution, where as introduction of denaturant such as GdmCl makes a polymer less 
compact. When plots of t^o vs rj/riQ are extrapolated to zero viscosity, the intercepts give the 
value of internal friction which increases with increasing value of kc and the slopes give the 
dependence of reconfiguration time on solvent viscosity. While compared with experiments 
nia, the results are consistent with experiments where the intercepts decreases as the 
denaturant concentration increases (low kc)- Typically the time scale due to internal friction 
calculated from the intercepts are in the range of 10 — 50 ns, which are also in good agreement 
with the experiments. Similar is the case with Zimm chain. 

Next we plot Tnm against \n — m\, the results obtained here are very fascinating. In case of 
CRIF, as shown in Fig. ([^, when the value of kc is small Tnm first increases with increasing 
|n — m|, until it reaches a hump, from where it starts declining. This happens due to two 
opposing factors. First one is MSD of the distance between the and monomer, for 
which an analytical expression is derived in the form of Eq. MSD for every pair of 

monomers at f = T^m is shown in Fig. (§. which clearly shows as the distance between two 
monomers become greater, the MSD is also higher which results in faster correlation loss of 
Rnm- If correlation is decaying fast Tnm will consequently become smaller. Another factor 
influencing reconfiguration is the distance between two monomers n and m along the chain. 
This can be seen from the normal mode description, Rnm is a sum of normal modes and as 
|n — m| decreases, the contribution of higher normal modes will be greater which relaxes 
much faster. Thus if \n — m\ is small Tnm should be small. Because of these two competing 
effects the reconfiguration time passes through the maxima [TB] . 

The trend of Tnm vs \n — m\ is quite different at higher degrees of compactness as shown 
in Fig. (|^. Again the plausible reason behind this might be the competing effect of the 
two factors, MSD and relaxation times of the normal modes describing the part of polymer 
between the and monomers. When the value of kc becomes higher and the polymer 
becomes more compact. If one looks very carefully into the expression for = 
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^R(Z) 


^cR(z)iF + %" 2 ip 2 ^ 2 , it seems with kc being independent of Qnt relaxation should be faster in 

presence of kc- But Cmt also becomes large with increasing kc and as told earlier ~ 

. 2 2 

is greater for higher modes in comparison to the lower ones, which is shown in 
Fig. ([^. In a nutshell slowing down of relaxation time is more for higher normal modes 
compared to the lower ones in presence of higher kc ■ So, at higher kc, compared to lower 
value, relaxation is more slower when the distance between and monomers is small, 
compared to the situation when distance between n and m is longer. MSD is reduced as 
well which can be seen in Fig. ([^ which is an obvious outcome of compact state of polymer. 
Because, in this new model beside the polymer being more compact, the movement of each 
monomer is very restricted. But it is clearly observed from Fig. 0 , the decrease of MSD is 
not consistent. With increasing distance between n and m, MSD is more controlled. That 
is why when |n —m| is small, effect of MSD become more prominent and Tnm decreases with 
increasing |n — m| and for higher \n — m\ the dynamics of the polymer between n and m is 
controlled by the the normal modes instead of MSD. As a consequence reconfiguration time 
increases with increasing \n — m\. 

When similar calculations are done for Zimm model, as can be seen in Fig. Q even 
for small kc reconfiguration time decreases with increasing distance between two monomers 
similar to the second trend discussed for CRIF with higher value of kc- The relaxation time 
of the normal modes describing Zimm chain are usually faster than that of Rouse [22] • That 
is why, at even smaller value of kc the second trend is be observed. This trend is very similar 
to the experimental results found by Nettels et al |16] where they estimated reconhguration 
time as a function of separating distance between two groups from FRET. If the value of kc 
is considered even smaller, the hrst trend observed for CRIF is recovered again which is not 
shown here. 

Looping time is estimated following the Wilemski Fixman (WF) approach, where a sink 
is placed between the two monomers which will eventually form a loop and has the closed 
form generalized expression given in Eq. (12). Fig. ([^ and Fig. (10) are showing the 
dependence of Tnm,loop on |n — m| for CRIF and CZIF respectively. Variations of Tnm,loop vs 
\n — m\ is presumably due to two competing factors as in the case of Tnm- 
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(a) 


(b) 



FIG. 1: Polymer model (a) Rouse, (b) RIF, (c) CRIF. 

V. CONCLUSION 

The work addresses the looping dynamics of single polymer within WF approximation 
including the effect of internal friction at various degrees of compactness. Looping and re- 
conhguration times are calculated using the recently proposed model compacted Rouse with 
internal friction (CRIF). Same set of calculations are carried out with compacted Zimm with 
internal friction (CZIF) where hydrodynamic interactions are accounted for. The strength 
of our method lies in its simplicity. We show that without carrying out expensive sim- 
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2 4 6 8 10 


P 


FIG. 2: Plot of 


vs p. See text for the values of the parameters used. 



p/il„ 


FIG. 3: Plot of end-to-end Reconhguration time vs rj/rio for GRIP. See text for the values of the 
parameters used. 
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FIG. 4: Plot of end-to-end Reconfiguration time vs r]/r]Q for CZIF. See text for the values of the 
parameters used. 



FIG. 5: Plot of Reconfiguration time vs |n —m| for GRIP. See text for the values of the parameters 
used. 
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FIG. 6: Plot of MSD at vs \n — m\ for CRIF. See text for the values of the parameters used. 



FIG. 7: Plot of Reconfiguration time vs |n —m| for GZIF. See text for the values of the parameters 
used. 
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FIG. 8: Plot of MSD at vs \n — m\ for CZIF. See text for the values of the parameters used. 



FIG. 9: Plot of Looping time vs |n — m\ for GRIP. See text for the values of the parameters used. 
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FIG. 10: Plot of Looping time vs |n — m| for CZIF. See text for the values of the parameters used. 

ulations even a simple model like ours confirm experimental results. It can account for 
different intercepts of reconfiguration time (r^m) vs solvent viscosity {rj) at different denatu- 
rant concentrations. Presently our model does not have a fluctuation-dissipation theorem for 
the internal friction m- Currently we are working towards a model having a fluctuation- 
dissipation theorem for the internal friction as well. Another aspect is the crowding on 
looping |2T]. It is expected as crowding increases internal friction should increase. This will 
be an interesting future problem to investigate. 
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